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Abstract We investigate the total kinetic powers (Lj ) 
and ages (iagc) of powerful jets of four FR II radio 
sources (Cygnus A, 3C 223, 3C 284, and 3C 219) by 
the detail comparison of the dynamical model of ex- 
panding cocoons with observed ones. It is found that 
these sources have quite large kinetic powers with the 
ratio of Lj to the Eddington luminosity (iEdd) resides 
in 0.02 < Lj/L-Edd < 10. Reflecting the large kinetic 
powers, we also find that the total energy stored in the 
cocoon (Ec) exceed the energy derived from the min- 
imum energy condition (£^min): 2 < E'c/^min < 160. 
This implies that a large amount of kinetic power is 
carried by invisible components such as thermal lep- 
tons (electron and positron) and/or protons. 
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1 Introduction 



individual (Cygnus A, 



The total kinetic powers of AGN jets is one of the most 
basic physical quantities of the jet. It is however dif- 
ficult to estimate Lj , since most of the observed emis- 
sions from AGN jets are of non-thermal electron ori- 
gin, and the electromagnetic signals from the thermal 
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and/or proton components is hard to detect. Hence, 
the free parameter describing the amount of the invis- 
ible plasma components always lurks in the estimates 
of Lj based on the non-thermal emissions. 



Recentlv. iKino and KawakatuI (|2005 1 (hereafter KK05)| 
proposed a new estimate of Lj for FRII radio galaxies 
based on the dynamical model of cocoon. Jets in power- 
ful radio galaxies are expected inflate a cocoon into the 
surrounding intra-cluster medi um (ICM) whi ch is over- 
pressured against the ICM (Be gelman and Ci ofE 1989^). 
From their model, Lj and iagc can be derived by com- 
paring the cocoon model with the actual morphologies 
of the cocoon based on the radio observations. How- 
ever, at the moment, this model has been only applied 
to Cygnus A. The extension of the number of samples 
are evidently crucial for exploring general characteris- 
tics of powerful AGN jets. For this purpose, we apply 
the method of KK05 (with a slight modification) to 
four FR II radio galaxies (Cygnus A, 3C 223, 3C 284, 
and 3C 219). 



2 Cocoon model 

Following KK05, we briefly summarize the dynamics 
of cocoon in over- pressured regime, namely Pc > Pa, 
where Pc and P^ are the cocoon pressure and the pres- 
sure of the ambient ICM, respectively. The basic equa- 
tions which describe the equation of motion along jet 
axis, equation of motion of the sideways expansion, and 
the energy equation are expressed, respectively, as 



— =/Oa(rh)VhW^h(i), 
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where wj, pa, Wh, Vc, Ah, and 7c are the velocity of 
jet, the density of ambient medium, the advance ve- 
locity of cocoon head, the velocity of sideways ex- 
pansion, the cross sectional area of cocoon head, and 
the specific heat ratio of plasma inside cocoon, re- 
spectively. In these equations we also assume that 
the jet is relativistic, and that Lj is constant in time. 
nM = Jl^^Vh{t')dt', rcW = Jl^^v,{t')dt\ and 
Vc{t) — |7rrc(t)^rh(t) are the distance from the jet apex 
to the hotspot, the radius of the cocoon body, and the 
volume of the cocoon where i,nin is the initial time of 
source evolution. The declining mass density of ICM 
Pa is given by pa(r) = Pa{r/n^)~" where fh = rh(tago)- 
A cartoon of the cocoon model is illustrated in Fig. [TJ 
We slightly improve the model of KK05 as follows; (i) 
the geometrical factor of T4 is more accurate, and (ii) 
the inclusion of PdV work done by cocoon. 
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Fig. 1. — A cartoon of interaction of the ICM with 
declining atmosphere and the relativistic jet. 

We further assume the lateral expansion of the co- 
coon as Ac{t) — Ac{t/tiigc)^ , where Ac{t) — 7rrc(i)^ is 
the cross sectional area of cocoon body. By assuming 
that all physical quantities also have a time-dependence 
of form A ~ A (i/^agc)^, Vc is determined as 



Vc{t) = Vc ( - — 
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From this, we obtain the fohowing expressions: 
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where Cvh = i(7c - 1)(0.5A:)-"[3 - (2 - 0.5a)X]/[4 - 
(1 - 0.5a)X], Cvc - 0.5X/^i/2, Cpc = (0.5X)", Cah = 
[X{-2 + 0.5a) + 3]"", and vq = fh/iago- Note that 
the only difference from KK05 is Cvh- Although it 
is not dealt in the present stu dy, the main differ- 



ence from the previous m odels (jBicknell et al.l Il997 



Kaiser and Alexandeiill997l ) is that these solutions can 
describe non-self-similar evolution. Note that when 
self-similar condition is imposed, time dependence in 
the physical quantities becomes identical to these pre- 
vious models. 

Throughout this paper, we assume 7c = 4/3 since co- 
coon is expected to b e dominated by relativistic parti- 
cles ( Kino et al. 2007t ). We also assume that the aspect 
ratio of the cocoon (TZ = rc/rh) to be constant in time 
(self-similar evolution). Since n{t) oc t^(2-5-o.5a)-3^ 
this assumption leads to X = 6/(5 — a). 



3 Comparison with the observation and the 
model 

By substituting the observed values of fh, fc (or 7?.), 
Ah, Pa, and a in Eq. ([5]) - (O, Lj, iagc, and also Pc are 
determined uniquely. In this section, we explain how 
to extract these quantities from the observations and 
show the obtained Lj and tagc 

In the left panels of Fig. [21 we show the VLA radio 
images of Cygnus A, 3C 223, 3C 284, and 3C 219 in 
logarithmic scale. Contours in linear scale (green lines) 
are also displayed for the purpose to determine the 
position of hotspot accurately. The overlaid straight 
lines which crosses perpendicular to each others at the 
hotspot are the lines we used to measure (black 
lines) and A-^ (red lines). Here, A-^ is measured as the 
cross sectional area of the radio lobe at the position of 
the hotspot. Contrary to A-^, it is difhcult to measure 
Tc from the VLA radio images since the cocoon emis- 
sion coming from the cross section perpendicular to the 
AGN core location is very dim at GHz radio frequency 
because of the synchrotron cooling. Therefore, we treat 
Vc = TZr^ as a free parameter. In order to take account 
of large ambiguity, we examine sufficiently wide range 
of 0.5 < 7?. < 1. In the present study, we neglect projec- 
tion effect on for simplicity. Under the assumption 
of self-similar evolution, Lj and tage depend on as 
Lj oc and iagc oc rjj, respectively. However, since 
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the angle to the line of sight is not expected to deviate 
largely from 7r/2 in radio galaxies, this effect does not 
cause significant change in our results. 

The adopt ed values of Pai ct, and Pa are based on the 



literatures oflRevnolds and Fabian ( 1996): ISmith et al 



tOO± for Cvgn us A. ICroston et all (l2004h for 3C 223 
and 3C 284, and iHardcastle and Worralll (|l999l ) for 3C 
219, respectively. Note that since most of the radio 
sources show asymmetries in their shape among the pair 
of the lobes (Fig. [2|) , we analyze each side of the lobe 
independently. However, we only analyze jet in the 
western side for 3C 219 since the eastern lobe showed 
severe deformation. 

The resultant Lj and tago are displayed in the right 
panels of Figs. [H The red lines displayed in the figure 
are the resultant range of Lj and iago, and their ranges 
reflect the uncertainty in TZ. The age and power depend 



oc 7^'*-" and Lj oc 7^ 



2c(-8 



on the aspect ratio TZ as 
and, therefore, satisfy Lj oc ^ago- Since a does not ex- 
ceed 4 in any of four sources, a lower aspect ratio corre- 
sponds to a higher power with a lower age. The range 
included in the shaded region is the forbidden region 
since the overpressure condition is violated. The Ed- 
dington luminosity (L-Edd) of each source, is also shown 
(green line) in these figure for compa rison. The black 



hole masses (A/bh) a re taken from |Tadhunter et al . 
12003) for Cygnus A, iMarchesini et all |2004h for 3C 



223 and 3C 219. For 3C284, we de rive Mr h from the 
B-band magnitude of host galaxy (jShi et al.l 120051 ) by 
using the empirical c orrelation of B-band m agnitude 
and black hole mass (jMarchesini et al.ll2004j ). In Ta- 
ble [Tl we summarize the allowed values of Lj and tage 
and the other relevant physical properties of the co- 
coon. Reflecting the asymmetries in their shape, the 
obtained Lj and tage show discrepancy among the pair 
of lobes especially in 3C 223 and 3C 284, and 3C 219. 
Since it seems natural to suppose that the properties 
of the jets are intrinsically symmetric and equal power 
and age on both side, we expect that the discrepancy 
is due to the asymmetry and/or inhomogcneity of ICM 
density profile. Here we interpret that the actual values 
of Lj and tagc is in the range obtained from both lobes. 



4 Discussions 

4.1 OnLj/iEdd 

In Tabledl the total kinetic power of the jet normalized 
by the corresponding Eddington luminosity, 2Lj/LEdd, 
is displayed. It can be seen that 2Lj/iEdd takes quite 
high value ranging from ^ 0.02 to ~ 10. Postulat- 
ing that the relativistic jet emanating from the AGN is 



powered by some part of released gravitational energy 
of the accreting matter (e.g., Marscher et al. 2002), 
these values give the minimum mass accretion rate nor- 
malized by Eddington mass accretion rate. Interest- 
ingly, our results indicate that mass accretion rates are 
super-Eddington ones in some FR lis (3C 219 and 3C 
284) since 2Lj/LEdd ^ 1- 

4.2 On the plasma content 

It is intriguing to explore how much amount of the 
internal energy, Ec — PcVi:/(7c — 1), is deposited in 
the cocoon compared with the widely-discussed mini- 
mum energy, i?min, of the radio lobe obtained from the 
minimum energy condition for the non-thermal elec- 
trons and magnetic fields. Here we calculate i^min 
base d on the observation of 178MHz band radio emis- 
sion I Hardcastle et al. 1998[ l. In Tabled] we summarize 
the resultant -Emin, Ec, and r/min which we define as the 
fraction of £^c to -Bmin (?ymin = Ec/Eynin)- In all sources 
Ec exceeds i^min, and the range of the ratio is obtained 
as 2 < rimin < 160. This implies that minimum energy 
condition is unlikely to be realized in these sources. 

Using 77inin) here we investigate the plasma content 
in the cocoon. Considering the components of energy, 
Ec is sum of the energy of the non-thermal leptons 
(electron/positron) and, if present, the energy of the 
unobservable particles such as thermal leptons and/or 
protons. The large excess energy of Ec from i^min is 
due to contributions from either of these components. 

Let us consider the possibility of large contributions 
from the non-thermal leptons. In this case rymin is given 
as rymin = /Umin whcrc U^'^ and J7min are the 
energy density the nonthermal leptons and the mini- 
mum energy density (C/min = -E-min/K;), respectively. 
It is useful to express f/^'^/f/min in terms of U^^ /Ub, 
where Ub is the energy density the magnetic field since 
U^'^ /Ub is widely investigated by a lot of authors (e.g., 
Kataoka and Stawarz; Croston et al. 2005). In the case 
of power-law distributed leptons, synchrotron luminos- 
ity (L^) is determined as cx U^^U^^^Vc- Hence 
with L^ and Vc observed, Uc and Ub are not indepen- 
dent of each other, and [/^"""/J/min can be written as 
U^^/Unun = 0.5{U^^/UBf^'^- Since recent studies 
shows that 1 < U^'^ /Ub < 10 on average, the value 
of 77min is expected to be up to 1 at most. Thus, non- 
thermal leptons are unlikely to be the main carrier of 
the energy. We conclude that significant amount of en- 
ergy is carried by invisible components such as thermal 
leptons (electron and positron) and/or protons. 
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Fig. 2. — left p anels: LoKarithmic-scal ed 5-GHz VLA map of Cyg nus A dPerley et al. 1984 ) and 1.5GHz VLA 
maps of 3C 223 (jLeahv and Perlevlll99l[ ). 3C 284 |Leahv et al.lll986l ). and 3C 219 |Clarke et aljfliii ) with linearly 
spaced contours (green lines) are displayed. The straight lines overlaid in each map denote the lines we have used 
to measure rh (black lines) and Ai^ (red lines), right panels: The red lines show the obtained range of Lj and iago- 
The shaded regions show forbidden ranges where the overpressure condition {Pc > Pa) is not satisfied. Also the 
Eddington luminosities (green lines) are displayed for comparison. 
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Table 1 The obtained properties of the jet and the cocoon together with minimum energy of the radio lobe. 



Source 
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21/j/LEdd 


E 


C 


E ■ 






(lO'*^ erj 


IS s ^) 


(Myr) 


(Mq) 






(10^0 


ergs) 


(lO^o ergs) 






Cygnus A E 


0.4 - 


2.6 


19 - 


47 


2.5x10^^ 


0.025 - 


0.16 


3.4- 


8.8 


1.4 


2.3 


- 6.1 


Cygnus A W 


0.35 - 


1.1 


30 - 


53 


2.5x10^ 


0.021 - 


0.068 


3.2 - 


5.7 


1.4 


2.2 


- 4.0 


3C 223 N 


0.15 - 


2.9 


140 - 


610 


3.2x10^ 


0.072 - 


1.43 


16- 


70 


0.88 


18 


- 79 


3C 223 S 


0.071 


- 0.2 


330 - 


560 


3.2x10^ 


0.034 - 


0.097 


6.9 - 


- 12 


0.88 


7.8 


- 13 


3C 284 E 


0.3 - 


18 


32 - 


260 


8.2x10* 


0.053 


- 3.4 


14 - 


110 


1.8 


7.7 


- 62 


3C 284 W 


0.1 - 


3.6 


100 - 


630 


8.2x10^ 


0.018 - 


0.67 


11 - 


68 


3.0 


3.7 


- 23 


3C 219 W 


2.6 - 


43 


37- 


150 


6.3x10^ 


0.65 ■ 


■ 10 


63 - 


250 


1.6 


40 - 


160 



5 Summary 

In this paper we have investigated the total kinetic 
power and the age of the rclativistic jet in four FRII 
radio galaxies (Cygnus A, 3C 223, 3C 284, and 3C 219). 
Below we summarize our main results. 

(I) A large fraction of Eddington power in the range 
of 0.02 — 10 is carried away as a kinetic power of the 
jets in the FR II sources. 

(II) The energy deposited in the cocoon, i?c, exceeds 
the minimum energy, i^mim by a factor of 2 — 160. 
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made use of SAOimage DS9, developed by Smithsonian 
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